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SUMMARY 

The spectrum concept is employed to study the selective effect of 
a stream contraction on the longitudinal and lateral turbulent velocity 
fluctuations of the stream. By a consideration of the effect of the 
stream contraction on a single plane sinusoidal disturbance wave, 
hiathematically not dissimilar to a triply-periodic disturbance treated 
by G. I. Taylor, the effect on the spectrum tensor of the turbulence 
and hence on the correlation tensor are determined. Lack of inter- 
ference between waves follows from the postulation of a low level of 
turbulence; this and the assumption of an inviscid fluid imply neglect 
of decay effects. The compressibility of the main stream is taken 
into account, but the density fluctuations associated with the tur- 
bulence is assumed to be negligible; this would be the case if the 
- turbulence originated from wakes and boundary layers in the very low 

speed portion of the flow. For an axi symmetric contraction and a 
particular isotropic initial turbulence some explicit results are 
obtained. The one-dimensional longitudinal spectrum is found to be 
distorted (as well as reduced in amplitude) with its peak shifted well 
to the right of the initial position above the zero of the wave-number 
scale. The selective effect of the contraction on the mean square 
longitudinal and lateral components of turbulent velocity is found to 
be given uniquely when the initial turbulence is isotropic, regardless 
of the details of the spectrum. If the initial Spectrum is anisotropic, 
as, for instance, that produced by a damping screen, then the selective 
effect is altered. 

In a crude extension, decay effects outside the scope of the 
theory are allowed for in first approximation. With this extension, 
a comparison with experiment is made of the selective effect on tur- 
bulent intensity where the estimated decay effects are comparable with 
the contraction effects. The agreement is good for the longitudinal 
component, very poor for the lateral component, the experimental data 
themselves being in conflict. 
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INTRODUCTION 

The generation of a wind-tunnel flow is always accompanied by a 
certain amount of turbulence; this is one respect in which the flow fails 
to simulate free-f light conditions. Measurements in the t unne l , particu- 
larly those sensitive to boundary-layer behavior, are known to be affected 
by this turbulence. Accordingly, the tunnel designer attempts to reduce 
'the intensity to the lowest practicable level. The use of honeycombs and 
damping screens in a large low-speed section (settling chamber) followed 
by a sharp contraction to the much-higher-speed working section is known 
to be effective. The honeycombs and screens located in a low-speed 
section reduce the absolute level of the turbulence with little drag 
penalty; then the relative level is greatly reduced by the large gain 
in tunnel speed through the contraction, aside from any effect of the 
contraction on the absolute level. 

Once the characteristics of honeycombs and screens are known, the 
further quantitative estimate of the reduction in turbulence involves a 
knowledge of the effect of the tunnel contraction 1 on the turbulence. 

It is known that the longitudinal component of the turbulence is greatly 
reduced (in absolute value) by the contraction; the behavior of the 
lateral component appears , on the other hand, to vary from no change to 
a substantial increase. Prandtl (reference l) obtained a quantitative 
estimate of the first effect by considering the conservation of energy 
for a perturbed longitudinal filament: if the initial stream speed 

is U, the filament speed U + u, with u<<U, and the final stream speed 
is ZU, then the final filament speed must be ZU + Z - lu; that is, the ' 
contraction reduces the longitudinal perturbation velocity u by the 
factor Z . For the lateral effect, Prandtl applied conservation of 
momentum to a small rotating cylinder of the fluid, with its axis 
cross stream, as the fluid traversed the tunnel contraction. He con- 
cluded that the lateral perturbation velocity v is increased by a 
factor -y/T. 1 

Prandtl' s considerations on the effect of a stream contraction 
were limited, as has been noted, to particular idealized "eddies". 

G. I. Taylor later (reference 2) attempted more realism by treating a 
mathematically defined model of turbulence which amounted to vortices 
in parallelepiped partitions arranged in a regular three-dimensional 


1 The considerations of this paper are not limited to a wind-tunnel 
contraction: they may be applied to any stream tube of varying cross 

section large compared with the scale of the turbulence. 
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array. The changes in vorticity on traversing the contraction were 
determined from a theorem based on conservation of circulation for an 
inviscid fluid; the corresponding altered turbulent velocity pattern 
was then calculated. The final result of the analysis consisted in 
expressions for the root mean square longitudinal and lateral turbulent 
velocity components u ' and v * downstream of the contraction expressed 
as. ratios of the corresponding values upstream. 

The initial condition of isotropic turbulence (mean values unaffected 
by rotation or reflection of axes) was approximated by specifying the 
vortex partitions to be cubical. For this case the reduction in the 
longitudinal component u‘ was found to vary more nearly like 1.5 Z _1 
than the value Z -1 suggested by Prandtl. No explicit result was 
found for. the variation of the- lateral component, however: the calcu- 

lations contained a free parameter. 

Taylor's results for the longitudinal component agreed fairly well 
with the experimental data then available, but ’ it is now considered that 
the measurements were made too close behind the screens for the screen- 
produced turbulence to have been isotropic. On theoretical grounds, the 
objection to Taylor's theory is threefold: first, the decay processes 

of turbulent mixing and viscous dissipation, which result in a reduction 
of the mean intensity with axial distance in the wind tunnel, are neg- 
lected; second, the assumed model of turbulence fails to exhibit the 
spatial and temporal randomness of actual turbulence; third, no choice 
of the parameters in Taylor's model corresponds to isotropy. In a sense 
all three objections apply likewise to Prandtl 's results: no model was 

employed in his considerations, and hence no distinctions between the ' 
effects of isotropy and anisotropy were made. 

The second and third objections can be removed by working, not with 
a model of turbulence, but instead with a Fourier integral representation 
of a random turbulent field. The integral represents a superposition of 
plane transverse sinusoidal waves of all wave lengths and with apparently 
random phases and planes of polarization. This aggregate of plane waves 
constitutes the (three-dimensional) spectrum of the turbulence. Only the 
statistical aspects of this spectrum will be known, not, for example, 
the detailed phase relationships. Mean square velocity components may 
be obtained by an integral of certain spectrum functions in which the 
phase relationships are suppressed; these functions are included in 
the "spectrum tensor" (reference 3). 

» Taylor's concepts may be applied to find the effect of a stream 
contraction on a single plane wave. The effect under the assumptions 
is linear; therefore the superposition implied by the Fourier integral 
may be employed to obtain the contraction effect on a field of turbu- 
lence. In particular, if the initial spectrum tensor is known the final 
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spectrum tensor is determined; the initial and final mean square velocity 
components then result from quadratures. Indeed, from the same informa- 
tion, changes due to the stream contraction in correlations of velocity 
at different points may he calculated: use is made of the fact that the 

correlation tensor is an inverse Fourier transform of the spectrum 
tensor (reference 3). 

Accordingly, the injection of the spectrum point of view into 
Taylor’s original concept of the contraction effect makes possible a 
more realistic calculation of the changes in mean square velocity com- 
ponents. In addition, it provides much more detailed information con- 
cerning changes in the statistical properties of the turbulence; that 
is, it provides the changes in the spectrum tensor and in the correla- 
tion tensor. 

The ideas just outlined are developed in the present paper. The 
first section is devoted to an account of turbulent spectrum analysis 
in. a form specially adapted to the analysis of the contraction effect. 

In this account, which is a generalization of a development in refer- 
ence 4, the role of the spectrum tensor is subordinated to that of the 
individual Fourier components (plane waves) in contradistinction to 
the customary treatment. This approach has perhaps an auxiliary merit 
in providing some better physical insight into the significance of the 
spectrum tensor, 

Next the effect of a stream contraction on a single plane wave is 
calculated by an application of Taylor's concepts. The treatment is 
slightly more general in that compressibility of the main stream is 
allowed for. The density fluctuations associated with the turbulence 
are assumed to be negligible; this would be the case if the turbulence 
originated entirely from boundary layers and wakes in the very low 
speed portion of the flow. Following Taylor, the problem is linearized 
by postulating a sufficiently weak turbulence so that the self-distortion 
of the turbulent eddies is small compared with the distortion imposed 
by the contraction of the main stream; this together with the assumption 
of an inviscid fluid implies neglect of the decay of the turbulence. 

In succeeding sections the spectrum and correlation tensors down- 
stream of the contraction are expressed in terms of the corresponding 
initial tensors. For the special case of an axisymmetric contraction 
and isotropic initial turbulence the ratios of the root mean square 
longitudinal and lateral velocity fluctuations downstream and upstream 
are obtained explicitly in terms of the parameters defining the con- 
traction. For a particular subcase where the initial isotropic spectrum 
tensor is specified, the corresponding 'one -dimensional* spectrums .(as 
would be recorded' by stationary hot-wire probes) upstream and downstream 
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of the contraction are calculated; the specification is such that the 
upstream, one-dimensional spectrum corresponds to experiment (refer- 
ence 5) in a number of cases of isotropic turbulence. 

Most of the calculated contraction effects are amenable to experi- 
mental checks either directly or indirectly. The available experimental 
data, however, are limited to the changes in the root mean square veloc- 
ities. A comparison with these experimental data is given with an esti- 
mated allowance for decay effects outside the scope of the theory. Design 
curves of the changes in the root mean square velocity components neg- 
lecting decay are included for engineering purposes. 


SPECTRUM ANALYSIS 


Representation of turbulence by superposition of plane sinusoidal 
waves . - Suppose q-j_, q 2 , q 3 . represent the components of velocity in a 

turbulent field; that is, q-^, q 2 , and q 3 vary in an apparently random 
manner in space and time, and the mean values *s q 2 = q 3 =0. Sub- 
ject to certain conditions, a snapshot of this field at any instant can 
be represented as a set of three-dimensional Fourier integrals 


qjxq, x 2 , x 3 ) = 



i (k 1 x 1 +k 2 x 2 +k 3 X3) 

Q a (ki, k 2 , k 3 )e dk^dk 2 dk3, 


where a = 1, 2, or 3 and the significance of k-^, k 2 , and k 3 will be 

brought out later. A continuous representation of the turbulent field 
is obtained by allowing the Q a to vary with time. 


It will be convenient to abbreviate the Fourier integral to 





(la) 


and to introduce the companion equation 



(lb) 


t 
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where k = k ] _ j) k 2 , k 3 ; x = x ± , x 2 , x 3 ; dx (x) = dx-jdx 2 dx 3 . The second 

equation allows, in principle at least, the coefficients Q a (k) to he 

calculated. Mathematically, q^x) and Q a (k) are termed three- 

dimensional Fourier transforms of each other; use will he made of this 
relationship later. 

The velocity components q^ are connected hy the condition of con- 
tinuity. In many cases of practical interest these turbulent velocities 
originate from boundary layers and the wakes of obstacles in flows of 
low subsonic speed, so that associated density fluctuations may be 
ignored ; this is still permissible when the turbulence so produced is 
transported by a high-speed stream. Thus the incompressible form of 
the continuity equation may be used and the result is 

% k l + ^2 k 2 + % k 3 “ 0 

This relation may be written more compactly as 

=0 ( 2 ) 
a 

Physical interpretation . - The amplitude components are com- 

plex in general. According to equation (lb), then, the requirement that 
the velocity components q a be real implies that Q a (-k) is the complex 

conjugate of O^k) „ If corresponding terms for k and -k in equa- 
tion (la) are paired their sum is thus equal to the real quantity 


2 (Re Q a ) cos (k*x) - 2(Jm Q a ) sin (k-x) (3) 

The imaginary parts -cancel in the pairing, which implies that they con- 
tribute nothing to the integral. Expression (3) represents a pair of 
plane standing waves, a cosine wave and a sine wave, with normals in 
the direction k = (k x , k 2 , k 3 ), where x = (x 1; x 2 , x 3 ) is the radius 

vector to any point. The vector k is termed the wave-number vector 
and its magnitude k simply the wave number; the corresponding wave, 
length is 2it divided by the wave number. Since k is perpendicular 
to the wave front, it is sometimes referred to herein as the ‘wave 
normal * . 

The continuity condition, equation (2), states that both the real 
part (Re Q a ) and the imaginary part (im Q a ) of the amplitude vector 
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Q = Q a = (Q-p Q 2 > Q3) ai ' e perpendicular to the wave normal k; that. is, 

both waves of expression (3) are transverse. For each wave any one of 
the parallel planes containing both the local velocity vector q and 
the wave-number vector k (which is perpendicular to Q and hence to 
q) is called the plane of polarization. The cosine wave (real part) 
and sine wave (imaginary part) may be polarized in different planes in 
general; the necessary and sufficient condition they be polarized in 
the same plane is that 


Re Re Qg Re 

Im Q-|_ ~ Im Qg ~ Im Q3 


( 4 ) 


Equations (lb) are now seen to represent a superposition of plane 
sinusoidal waves (Fourier components) with all orientations of the wave- 
fronts (all directions of the wave-normal k) and all wave lengths (all 
wave numbers k) . Each wave is transverse, and all planes of polariza- 
tion are permitted. For each value of k there exists a cosine wave 
and a sine wave; their respective amplitudes and planes of polarization 
are different in general. The complex amplitude components O^k) 

express, in their real and imaginary parts, how the respective ampli- 
tudes and planes of polarization vary with the wave-front orientation 
and the wave number. 

Mean values; the correlation tensor . - Consider the spatial"*" 
mean value of the product of the velocity component q a at x and the 

velocity component qp’ at x' = x + r as x varies but the separation 

r of the two points remains fixed during the averaging process; this 
mean value is called a velocity correlation and is given the symbol 
R a p(r). There are nine such correlations, corresponding to a = 1, 2, 3, 

P = 1, 2, 3. The form R a p(r) has been shown to transform like a 

second-order tensor and has been designated (sometimes divided by q^) 
as the "correlation tensor" (reference 6). 


-•-If the statistical properties of the turbulence are independent of - 
position (homogeneous turbulence) and time-independent, an average at a 
given time over all space equals- an average at a given point (or pair of 
correlated points) over all time; a proof is given in appendix B. If the 
statistical properties vary slowly with time the space average will still 
approximate a time average over an interval just long enough to smooth 
out the fluctuations. 
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form r a p(k) defined by Batchelor (reference 3). as the Fourier transform 
of RopCr). The form r a p(k) is known as the spectrum tensor. The 

Fourier transform relations connecting the spectrum tensor and the corre- 
lation tensor are summarized as 

r o0(i£) = 


By use of the Fourier transform relationship, Batchelor demonstrated 
that is a second-order tensor and obtained a number of its proper- 
ties. Thus, for example, T a p is complex, in general, with = 1^*^ 

and the diagonal elements are real; also, r a p(-k) m 

It is of interest to observe that these same properties result immedi- 
ately from the identification of Q a Qp*. with Thus 

is complex, in general; Q^Qq* equals [Q a Qp*] *; and Q a Q a * is, 
of course, real. Furthermore, since Q^-k) = Q a *(k), Q a (-k)Qp*(-k) 
equals Q a *(k)Qp(k); hence ^(-k) = Tp a (k) . 

The distinction between cases where T a p is real and cases where 

it is complex may be given a physical interpretation. The product 
Q a Qp* > and correspondingly r a p, is seen to be real when the condition 

equation (4) is satisfied. This implies that the cosine wave and sine 
wave associated with wave number k are polarized in the same plane.' 

The alternative condition where QgQp*, and hence are complex 

implies polarization of corresponding cosine and sine waves in different 
planes. The velocity pattern of such a pair of waves is quite inter- 
esting: successive velocity vectors along a line in the direction of 

the wave normal k turn progressively about this line in spiral fashion; 
the tips of the vectors trace out a helical curve on a cylinder of oval 
cross section. 

Energy spectral density . - Each of the diagonal elements r lx , 

F 22 > ^33 °f the spectrum tensor may be interpreted as an 

energy spectral density. Thus, according to equation (6b) 
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B-q(O) = U]_ =/ / / r l:L (k)clT(k); therefore 


the differential 

energy component 
within the range 


|r n dt (k) represents the contribution to the kinetic 

1 per unit mass made by waves with wave number 

2 

dx(k). 


One-dimensional spectrum . - The elements of the three-dimensional 
spectrum tensor are not directly measurable; they may be. obtained by 
taking the Fourier transform of the measured correlation tensor. A hot- 
wire probe placed in the moving stream will, however, develop a fluctu- 
ating output voltage whose (one-dimensional) frequency spectrum (refer- 
ence 7) is related to a diagonal element of the three-dimensional spec- 
trum tensor. Thus by equation (6b) the contribution to the mean square 


velocity component q_ 2 (=E aa (0) ) from all waves with wave-number com- - 
ponents in the direction of the x-^-axis between jk^j and |kj_| + JdkJ 
is «o 


F a (k 1 )dk 1 


2 {IP 


r aa! ik 2 dk 3 J dk l 


(7) 


the factor of 2 accounting for suppression of negative values of kj_. 
■The function F a (k 1 ) is the one-dimensional spectrum corresponding to 
the velocity component q. a ; the values a = 1, 2, 3 correspond respec- 
tively to the longitudinal and two lateral spectrums. The particular 
spectrum obtained depends on the arrangement of the hot-wire probe 
elements. 


EFFECT OF STREAM C0NTEA.CTI0N 

\ 

Consider now that the turbulent velocity pattern q.^, q.g, q.^ is 
carried along by an inviscid general stream with velocity U(x^) in the 
x-^-direction. Consider also that q.^, q^, q. 3 are so small that their 

effect on the streamlines may be neglected as the flow traverses a wind- 
tunnel contraction. The contraction will, however, distort the shape 
of fluid elements. (See fig. 1.) The vorticity distribution will be 
forced to alter accordingly to conserve the circulation about each 
element. The net result will be an altered pattern of turbulence. Each 
plane wave (Fourier component) Q-j_, Qg.> Q3© — will, in fact, be 
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altered independently under the linearizing assumption to "be made; the 
over-all effect on q-^, <3.2, q.3 can be obtained by the summation expressed 

by equation (l). It thus suffices to consider the effect of the contrac- 
tion on a single representative plane vave„ 


Effect of Contraction on Representative Plane Wave. 

Velocity and vorticity at upstream station . - Designate by A a 
reference station upstream of the contraction and by B a reference 
station downstream of the contraction. (See fig. 1.) Let a typical 
Fourier component (plane wave) of the turbulent field q a (a = 1, 2, 3) 

at station A be represented at time t = 0 by 

C = Ce 1( ^ (8) 

This rave, equation (8), is supposed to he carried along h j the main 
stream, with velocity U 0 

The vorticity w a is obtained from the curl of equation (8) as 



0,7 


where 6 ap^ 


0, if any pair of subscripts are equal 
< 1, if afiy are in cyclic order 
-1, if c $7 are in anticyclic order 


( 9 ) 


Distortion of a fluid element in passage through contraction . - 
Suppose the contraction is such that the stream velocity U is increased 
by a factor between stations A and B while the breadth and 

height of the tunnel are reduced by factors %2 and respectively. 

(See fig. 1 (a).) In traveling from A to B an initially cubical 
element of fluid of edge D will be distorted into a parallelepiped of 
edges Z-|D, Z2D, Z3D (see fig. 1(b)); a particle in the element originally 

(t = 0) a vector distance x from a comer particle will finally (t = t) 
be found a distance from the comer particle, where £ is related 

to x by ~ 
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Ci = l]*! 

C 2 = z 2 x 2 > 

^3 = Z 3 X 3 


( 10 ) 


In this argument the modification of the streamlines due to the turbulent 
velocity fluctuations has been neglected. This implies that the relative 
displacement of two adjacent particles due to the superposed turbulent 
motion is small compared with the displacement due to the tunnel contrac- 
tion. This key assumption, due to Taylor (reference 2), linearizes and 
vastly simplifies the problem. The limitations imposed by the assumption 
are discussed later under "Decay Considerations". 


The velocity ratio 2 and the lateral and vertical contraction 
ratios and l 3 are related by the continuity condition 


a2i2 2 2 3 =1 

where cr is the ratio of stream densities at stations B and A; the 
density is considered uniform at each station in accordance with the 
initial assumption of negligible turbulent density fluctuations. 

Yorticity at downstream station . - The vorticity is carried along 
by the flow, the fluid elements undergoing the distortion pictured in 
figure 1(b), to the approximation used. During the motion the strength 
changes in such a way as to maintain the constancy of circulation of the 
fluid elements. The changes are expressed by the equations for the 

transport of vorticity in the Lagrangian form, due to Cauchy (see 
reference 8), 

s„ B . i 

a 


E- 


m A 

P Sx R 


-*-This statement is exact for the postulated inviscid fluid. The 
modification produced by the diffusive effect of viscosity, in the case 
of a gas, becomes appreciable for the smaller eddies or higher wave 
numbers; for this analysis a criterion for neglect of viscous effects is 

vk^<<|dU/dx|. (See "Decay Considerations", equation (43).) 
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where a is the density ratio between stations B and A, and the deriv- 
atives /dx ft express the effect of the fluid distortion. 1 Evaluation 

GO p 

by means of the distortion equations (10) yields simply 


or, in expanded form 





ol^ 

al^ 2 A 


> 




(11a) 


These equations relating downstream and upstream vortioity embody the 
entire dynamics of the contraction effect. The equations are not limited 
to the plane sinusoidal waves discussed earlier, but apply to any (weak) 
vorticity distribution whatsoever. 


The above derivation of the vorticity changes is substantially in 
the form given originally by G. I. Taylor (reference 2) for the case 
a = 1 (incompressible flow). In order to assess the influence of the 
simplifying assumptions a more general derivation based on the Navier- 
Stokes equations is given in the section entitled "Decay Considerations". 

By virtue of equation (9) as applied to (lla) the vorticity at 
station B is obtained explicitly as 




~ B 
0) 

a 


= lCTl a / /apT^py 1 -*- 
P,7 


(lib) 


where, it -will be remembered , x. is the radius vector to a fluid particle 
at time t = 0 when the fluid element is at station A in the moving 
coordinate system of figure l(b')« The corresponding vector to the 
particle at time t = t, when the fluid element is at station B, is £ 
in that figure. When equations (10) are used to express x in terms 
of 4 , the exponential h-x becomes , in expanded form, 


^These equations refer to axes moving with some fluid particle 
rather than axes fixed as in reference 8; the form of the equations is 
unaffected . 
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k*x 


k l^l k 2^2 k 3^3 

J 1 + l 2 J 3 


The right-hand side may he expressed as where 



_ ^1 ^2 
= V V 


^3 

l 3 


defines a new wave-number vector. 


Telocity at downstream station, general case . - The velocity dis- 
tribution whose curl in the £ 2 , S3 system is given by equation (lib) 

and which satisfies continuity is found to be expressible in the form 


with 








VVa 

7 2 „ 2 
2p x j 


( 12 ) 


(13) 


where x is the magnitude of the wave-number vector x. This result 
is the general solution for the contraction effect on a single plane 
wave. 


Equations (12) and (13) admit of a simple (but not obvious) 
geometrical interpretation: traversal of the stream contraction alters 

the initial plane wave, equation (8), so that its wave-number vector 


k = (kq, k 2 , k 3 .) is transformed into x = ~ , J~, — and its amplitude 

pa R « T K> * 2 

vector (Qq , Qg , % ) is transformed into the projection of 

(Qq A Aq, §2 A /*2' ^5 A /^ 3^ on a plane normal to the new wave-number 
rector x. 


Telocity at downstream station, axisymmetric contraction . - In 
case the stream contraction is axisymmetric 1 a considerable simplifica- 
tion results. The condition for axisyrametry 2 2 * 2 3 , with use of the 

continuity equation (2.2), reduces equations (13) to 


A contraction such that all cross sections of the tunnel are 
similar, whence 2 2 = Z3 > is termed axi symmetr i c ; the sections need not 

be circular. 
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§2 B 


% B 


where e = Z 2 /Z 2 . 


~ A , 2 ,2 ,2 

Ql + k 2 + k 3 

*1 ek-i 2 + ko 2 + k. 


_1_ 

l 2 


_1_ 

*2 


x A Qik^Ci-e) 

2 + ek 2 + k 2 + k . 


Q^kgd-e) 

“i 2 + k 2 2 + *3* 




0 

] 


J 


(14) 


The considerably greater complexity of equation (13) is perhaps obscured 
by the purposely expanded form of equation (14). 

If the initial wave normal k is perpendicular to the (longitudinal) 
X]_-axis , the component k-]_ vanishes and equations (14) reduce further to 


% B = % A Ai 

q 2 b = q 2 a/z 2 


(15) 


q 3 b = Q 3 A A 2 

** A 

The same equations result when may be neglected in comparison with 

Q 2 A and Q 3 ^ , ' that is, when the amplitude vector is substantially normal 
to the x^-axis. Equations (15) state that an axisymmetric contraction 


defined by l Z 2 alters these waves by a factor of 

in the longitudinal velocity component and a factor of -s— in the lateral 

t _2 

velocity components. These equations apply only to particular types of 
waves; yet when the contraction effect is later integrated over the random 
aggregation of waves representing isotropic turbulence the over-all 
results are found not to differ greatly from the simple factors 

and respectively. 

h. ]2 

The same factors were obtained by Prandtl (reference l) for other 
special disturbances: the factor -i- from energy considerations for 
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purely longitudinal disturbance velocities, and the factor 


from momen- 


tum considerations for a rotating cylindrical element of fluid with axis 
normal to the stream. 


Effect of Contraction on Spectrum and Correlation Tensors 

Effect of contraction on correlation tensor . - The analysis herein 
lead s first to the changes in' the spectrum tensor r a p A (k) •* r ap B (x ) . 

Then the corresponding changes in the correlation tensor may be 
obtained from the Fourier transform relationship, equation (6b): 


CD 


(16a) 


E a^ B (2l) 


y =ffS r ' 




ix*r. 


(16b) 


In succeeding paragraphs r a p B (x) will be determined in terms of the 
initial spectrum tensor T^Ck) for various cases. 

Spectrum tensors at upstream and downstream stations in terms of 
the Qq . - In an earlier discussion the Fourier coefficients were'' 

chosen so as to define a field of turbulence confined to a large 
parallelepiped of volume t } and vanishing everywhere outside; for this 

JQ. 

^-_>oo -j- QoQb* was Be identified with the correlation tensor 


case 

r a|3* For station A. upstream of the contraction it will be convenient 

to specialize this parallelepiped to a cube of edge D. Such a cube 
will, however, be distorted into a parallelepiped of edges 7-jD, 2 2 D, 

ZjD by the stream contraction by the time it reaches station B down- 
stream. (See fig. 1(b).) The spectrum tensors for stations A and 
B, respectively, are therefore 


r op A te> 


Lim a, ’ b*. . 

IT*- 00 d 3 (^)^a ( F ) 


r o B 
L ap 


(x) = 


_ Lim Sjf 5 „ B/ ^ B*/ ,, 

(l) 


.(17) 
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Evaluation of spectrum tensor at downstream station, general case , - 
The identifications made in the last paragraph allow the spectrum tensor 
to he evaluated at station B in terms of the spectrum tensor at 
station A and the parameters defining the stream contrac- 
tion between stations A and B. For a single plane wave L = Q^e*— , 

which is transformed by the contraction into q^® = Q^e*— £, equa- 
tions (13) give 


** "R ~\ 

-T- 


►a 


2 A 


w-'E 




In the Fourier integral q^" is to be interpreted as dq^, q^® as 
dq^®, Qq^ as Q a ^dT(k) , and Q a ® as Q^dT (x) . Accordingly 


q^B = 2 1*2*3 


X 


a 




a. y Qp% k a 


, 2 2 

x 


(18) 


since = ^ T ( k )A*- T ( x ) • Thus 


t 2-j 2 

0 b q b* _ *1 *2 , }?_ 

Qa *a*p 


n An A* \ /" Qc^7 A *kykp Qp^Qy^ka Qy^^kykskakp | 

!“ p ’ir V v 2 * 2 ' y 2 * 2 * yV*‘ / 


The corresponding relation between the post-contraction and. precontrac- 
tion spectrum tensors is, by virtue of equations (17), 

(19) 


*1*273 


' * a*p 




t \ /ra7 A (k)krItp r 7P A (P) k r k a. r 76 A (k)k r liskaJ£p ) 

_ r< * yV^ 7 j 


where k is related to x "by 

k l> k 2> k 3 = l l x l> J 2 X 2> *3 X 3 ( 20 ) 

Special case : axlsymmetric contraction but arbitrary Initial 

spectrum . - When the contraction is axisymmetric (7^ * 2~E7 tbe squat i on 

of continuity in the form 
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E v> A * ° 

7 


( 21 ) 


may be used to simplify equation (19). The result may he written 


op ^ 


(x) = 


Z 


j/2 


aj3 ^ 


(k) + 


fraA 


P + ^lpAi-VlC 1 " 6 ) r 6 ) 2 

2 ! “b- _ 2 _ , _ 2 / p P p\2 

2 + k 2 2 + k 3 2 ) _ 


ek-^ + k 2 + k 3 


K* 


( 22 ) 


2 / 2 

where the ratio e = ^2 /h > f or a large speed gain in the contraction 
e « 1 . 

Special case: axisymmetric contraction, isotropic initial spectrum . 

A further simplification occurs when the turbulence at station A is 
isotropic. In that case, the spectrum, tensor r a o B (x) downstream of 
the contraction can be expressed explicitly to within an unknown multi- 
plicative factor G(k). This results from the fact (reference 3) that 
r o0 A (£) then be an isotropic second-order tensor; the isotropic 

property together with the continuity condition, equation (21), requires 
Tap (k) to be of the form 


r ap A (^) = G(k)(k 2 5 a ^ - k a k p ) (23) 

where 


>ap 


-G 


for a = p 
for a ^ p 


The right-hand side of equation (22) may be evaluated by means of 
equation (23). The diagonal terms reduce to relatively simple forms’: 


r il B <i) 


w 

7 2 
l l 


G(k) 


f* 2 - n 2 ) 


(eki + k 2 + fe 


1 


(24) 


r 22 ^ 


•lil 


1 6 2 


G(k) 


k^ - k 9 


2ki 2 k 2 2 (l-e) 

, 2 ~Z ~2 
eki + k 2 + k 3 


*!% 2 (k 2 - kx fi )(l-«) 2 


ek-, + k? + ki 


f j 


(25) 
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B B 

and r 33 "(x) is obtained from. F gg (*) ty replacing kg ty k 3 and 

vice versa. The relation, equation (20), between k and x applies 
here . ' 

One -dimensional longitudinal spectrums . - If the form of the initial 
spectrum tensor r op^(k) is known, the corresponding one-dimensional 
spectrums F^Ck^) can be calculated, according to the defining equa- 
tion (7) as applied at the upstream station A: 



kg, k 3 )dkgdk 3 


(26) 


A particular case of isotropic turbulence is of special interest 
(reference^): in equation (23) for T a p A (k) the function G(k) 

taken to be N(k 2 +'y2)“3, where N , y are constants. Then 


is 


•if - 


k 


ii fk 2 2 1 k 3 2 ) 


+ kc> + k-z 




and after integration 



This one -dimensional longitudinal spectrum is of the same fom as an 
empirical relation obtained in reference 5 for that of isotropic 
turbulence in the initial period; this agreement is the special virtue 
of the form assumed for G(k) . 


The one-dimensional lateral spectrum functions corresponding to 
the same G(k) are readily evaluated; they are 


F 2 A (k 1 ) = F 3 A (k x ) = 


cN^k-, 2 + y 2 ) 

; (*i 2 + r 2 ) 2 


The equality of the Fg and the F 3 functions results, of course, from 
the isotropy of the turbulence. 


The effect of the stream contraction on these one-dimensional 
spectrums is found by employing the post-contraction value of f , 
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B 

that is, the value appropriate to the downstream, station B. 

Since is a function of the local wave-number vector x at 

station B, the equation corresponding to (26) is 

o» 

** “ 2 XT I '“‘ B( V V x 3> dx 2 ax 3 

CD 

which is a function of 


For performing the integration and making later comparisons of 
spectrums it is convenient to transform from x-^, Xg, Xj to k^, kg, • 

kj, where 


X 1 = k l/*l 
x 2 = k 2/ Z 2 
*3 = k 3/ Z 3 


and to define F a ^( k i) = F a ®, such that 

thus _ 

F, 


^ *<? (*!>«*! -f 


• iXhJJ Ta * (*i' h’ h) 


dkgdk 3 


(27) 


B ' 

The spectrum tensor elements following an axi symmetric contraction 

have been evaluated in equations (24) and (25). With these values inserted 
and G-(k) specified as before, the integrations of equation (25) are best 
effected in polar coordinates. The results are expressed most simply in 
terms of a "normalized” longitudinal wave number k as incorporated 


in the two parameters 



t 



1 


The final result for the one -dimensional longitudinal spectrum following 
an axisymmetric contraction (lg = 2j) is 
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»i*( k i> = -rfs 3 + 4t + 12 + 

7 * _ 


+ 4s + 2t + st + 


The corresponding result for the one-dimensional lateral spectrums 
following an axi symmetric contraction is 


V = F 3 - 


JtN / (3s-2)t^ 


27 2 2 2 2 t 2 l s 2 


- 4(1-6) ( 


J 2s+t s+t / e A 1~6 f" 6s+5t (s+t) (3s+t) ( a 

+ T" l0g e (jTtJ + + t lo Se (iftjjj 


and for e«l (large speed gain) a simple "but very close approximation 




1 + 2k x 2 /7 2 


2 h Z 7 2 (l + k 1 2 / 7 2j 


+ 0(e, e In e) 


(The corresponding approximation for F^Os^) is not simple enough to 

warrant noting.) The parameters l 1} Z 2 = ^3 and c = I? 2 /Z} 2 are 

related to the initial and final Mach numbers of the main stream by the 
equations 


1%/ 


2 / 5 + M a 2 

l s + Mb 2 


2 m a { 5 + Mb 


% \5 


ma\ 3 (5 + Mb 2 


5 


For incompressible flow (Mg, M A -♦ 0) these reduce to 


ll = V U A 

2 = Z r1 / 2 

*2 L 1 


e = Z- 


(30a) 


(30b) 
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These post -contraction ape ctr urns, equations (28) and (29), are com- 
pared with the initial spectrums in figures 2 and 3, respectively; the 
comparison is "based on an assumed initial stream Mach number of 0,05 
(station A) followed by an axisymmetric contraction such that the final 
Mach number is 2.0 (station B); the corresponding parameters are 
l 1 = 29.8, = 0.382; e = 0.00016. Consider first the longitudinal 

spectrums, figure 2. Normalizing factors are incorporated such that 
the areas under the two curves, if replotted on a linear scale, would 
be the same; this normalization serves to differentiate changes in 
shape from changes in amplitude. The figure exhibits a rather striking 
distortion of the spectrum after traversing the stream contraction: 
the peak spectral density is shifted from zero wave number to k - 1.4 

along with a general shift of density to the higher wave numbers. Asso- 
ciated with this chang e in shape is a reduction in amplitude by the 

factor u b 2 /u a 2 , u a 2 and u B 2 being the respective integrals of the 

spectral density curves. These integrals are evaluated in a later 
section. 


The corresponding comparison for the lateral one-dimensional 
spectrums is made in figure 3. In this case the axisymmetric contrac- 
tion has made very little distortion in the spectrum. There is again 
a change in magnitude (this time an increase) in the ratio 




The changes in magnitude (that is, the changes in area under the 
spectra l density curves) correspond to the changes u B 2 /u A 2 and 
v b 2 / v a 2 in the mean square components of turbulence and are, at least 

qualitatively, well known. The predicted changes in the shape of the 
spectrum curves are apparently new. 


In the above comparisons both pre-contraction and post -contraction 
spectrums have been expressed in terms of the pre-contraction longi- 
tudinal wave number k l> whereas the local post-contraction wave number 

is = k i/li* Consider, however, a representative longitudinal wave 
which has the form cos k^x at station A and cos x-j_£ at station B. 

If x and £ are identified with the respective distances swept in 
time t by the moving waves over stationary hot-wire probes at 
stations A and B respectively, then k^_x = k^I^t and x-^£ « x^IJgt = 

ki 

— Z x U A t. Thus the (temporal) frequency seen by the hot wire in both 

cases is k x U A /2ito The comparison based on k x therefore constitutes, 
in effect, a comparison of the time spectrums that would be seen by 
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stationary hot-wire probes, in contradistinction to the space spectruxns 
discussed in the earlier parts of the paper. 


Effect of Contraction on Mean Square Telocity Components 
for Isotropic Turbulence 


The mean square velocity components of the turbulent field may be 
identified as the diagonal terms of the correlation tensor R a p (r ) with 
r set equal to zero. Thus 

u 2 = H u (0) 


v 2 = R 2 2(°) 
w 2 = 1*33(0) 


where u, v, w have been written for q-^ q 2 , q 3 , respectively. The 

evaluation of these means is much less laborious than the evaluation 
of the general correlati on tenso r. In particular, the evaluation of 

the ratio of the means Ug 2 /u A 2 , etc., may be made when the initial 

turbulence is specified to be isotropic but no further details of its 
spectrum are known. These ratios will be calculated in the following 
paragraphs . . 

p p 

Evaluations of u and at upstream station ,, - According to 

equation (6) 


R n A(0 ) = u a 2 



11 A (k)dx(k) 


For isotropic turbulence T a p has the form specified in equation (23), 
whence 

00 

U A 2 = J'g( k)(k 2 - k 1 2 )dv(k) 

— oo 

where G(k) is an arbitrary function. It is convenient to transform 
to spherical polar coordinates: 
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Then 


= k cos 6 


= k sin 6 cos <p 


= k sin 6 sin <p 


dT(k) = 


k^ sin 6 d 6 d<p dk 


u A 2 = 


k 4 G(k)dk d cp sin^e d 9 


For the present purpose the function G(k ) , which, together with the 
condition of isotropy, defines the turbulence, may he left unspecified; 
the integral involving G(k) will cancel out in forming the ratio 

u^/u A 2 . Let this integral have the value H; then 

u a 2 = | JtH 


By virtue of the assumed isotropy 


V A 2 = v k = f rtH 

O* O* 

Evaluation of ratio of u at downstream station to u at 

* O 

upstream station . - The mean value ujj is obtained from an integration 
involving the spectrum tensor after the latter has been, transformed by 
passage of the flow through the tunnel contraction; according to 
equation (16b) 


Eii B (0) 




r il B (x)<Mx) 


For the present case, where the spectrum tensor at station A is 
assumed isotropic and the contraction is axisymmetric, . the transformed 
tensor T]_}. S (x) has "been determined in equation. (24) 0 Thus 


7T_ W 

U B ~ , 2 


* G(k)k 4 (k 2 - kjjdx(x) 

( ek l 2 + k 2 2 + k 3 2 ) 2 
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Because of the unspecified function G(k) it is convenient to change 
the variables of integration from the components of x to "the components 
of k. In other words, a transformation is made from the "wave -number 
space" of station B to the "wave-number space" of station A. The trans- 
formation follows from the Cartesian relations 


dx (k) = dk-^dkgdkj 

dx(x) = dx 1 dx 2 dx 3 

dkj dkg dk-j 
1 1 ^2 ^3 


(33) 



Again the polar-coordinate transformation (equation (3l)) is made, with 
the result 


U B 


• 2n Jt 

= ^ k 4 C(k) dk ^ 


sin% d© 


Jq (e cos 2 © + sin 2 ©) 


•The first two integrals occur also in u A 2 (equation (32)), and they 
cancel in obtaining the ratio Uj^/u^ 2 ; thus 

2 


^B 

U A 


sin°0 d© 


2 ^1. <Jq (e cos 2 © + sin 2 ©)^ 


The final result may be written 


u- 


L B 


U A 


2 42 n 


-1 

1-6 


2-e 


7 ^ 


( 1 - 6 ) 


tanh”^- -y/l-e 


(34) 


and an asymptotic expansion for small e is 
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1 + |- e + (1 + e)ln ^ + 0(e 2 In e)j 

Equation (34) gives the ratio of the mean square longitudinal 
velocity fluctuation downstream of an axisymmetric tunnel contraction 
to the corresponding mean square upstream of the contraction, when the 
initial turbulence is isotropic. The contraction is characterized by 
an increase in the stream speed in the specified ratio 2^ and a 

decrease in the lateral dimensions in the specified ratio Z 2 ; tlie 
parameters 7 Z 2 , and e = Z 2 2 /Zq_ 2 are completely defined by the 
initial and final Mach numbers of the stream according to equations (30). 


The variation of u^/u^ with the speed ratio is plotted 

in figure 4 for two examples; in the first the flow is assumed com- 
pressible with a Mach' number 0.05 at the start of the contraction; in 
the second the flow is assumed incompressible (m a , Mg-»0). The Mach 

number scale at the bottom applies only to the compressible case, the 
?! scale to both cases. The salient characteristic of the curve is 

the marked reduction in the longitudinal component of turbulence with 
increasing speed ratio 2]_. 

Compressibility is seen to have but a secondary effect, which is 
appreciable only at supersonic speeds. Note (equations (30) and (34)) 
that with l 1 as the independent variable, the effect of compressibility 

appears only in the parameter e. The physical significance of 6 
follows from the definition of as the speed ratio provided by the 

p 

contraction and Z 2 as the area ratio of the contraction (in the 
axisymmetric case considered), with e = Zg 2 /^ 2 . Eor supersonic final 

speeds it is more proper to speak of a converging-diverging .nozzle than 
a contraction, the term "contraction" having been retained herein 
primarily for reasons of past usage. ' 

The basis of the compressibility effect may be summed up in the 
following way. The influence of an axisymmetric stream contraction 
arises from distortion of the fluid elements, as described by the 
parameters l 1 and Z 2 . (See fig. 1(b).) These parameters are 

related by the continuity condition oZ^Z 2 2 = 1,- where o is the 

density ratio. Thus compressibility, in allowing a to deviate from, 
unity, changes the relation between Z]_ and Z 2 somewhat, and con- 
sequently modifies the contraction effect. 
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The graph of equation (34) in figure 4 is primarily for illustrative 
purposes; a form more useful for engineering applications is given in 

figure 5. The single curve provides the variation of /\] ug^/u^ with 
both 2 ^ and e; Z]_ and e may he determined from ’the initial and 
final Mach numbers by means of the simple relation (30a) . 

Evaluation of ratio of v^ at downstream station to v^ at 
upstream station . - The value of Vg2 results from an integration 
involving the transformed spectrum tensor, according to equation (16b) 

OD 

*22 B (°) - V B 2 =jjj' r 22 B & dT & 




For isotropic initial turbulence and an axisyrmnetric contraction the 
transformed spectrum tensor iias evaluated in equation (25). 

Thus 



2k l 2 k 2 2 ( 1 “e) 

ek-j _ 2 + k 2 2 + k 3 2 


+ 


k 1 2 k 2 2 (k 2 - k 1 2 )(l-e) 2 
(ek/ + k 2 2 + k 3 2 ) 2 


dx(x) 


Again it is convenient to transform from x -3 P ace to k-space (equa- 
tions (33)) and to introduce polar coordinates k, <p, and Q (equa- 
tions (3l)). The integrations with respect to k and tp are readily 
disposed of, with the result 


H 




sin 3 e de - 2it(l-e) 



sln 5 e cos^e d6 
sin 2 0 + e cos 2 0 


+ 



alnSg cob 2 6 d8 
sln 2 0 + e cos 2 0 


where H = 


k^G-(k)dk, as before. Upon carrying out the integration 


2 8 


and dividing by v^_ = — itH there is obtained finally 


V B 


2 8Z 


2 


2 L 


2-6 e^ j_ , -1 r r * " 

77 - 7777 tanh -v/ 1-e 

1 “ e ( 1 - e ) 3 / 2 ^ 


( 35 ) 
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For small e this has the asymptotic expansion, 




+ e 


+ i 



Equation (35) gives the ratio of the mean square lateral velocity 
fluctuation downstream of an axisymmetric tunnel contraction to the 
corresponding mean square upstream of the contraction, where the 
turbulence has been assumed to be isotropic.- The variation of 

\/ v B 2 / v A 2 the speed ratio 2i is plotted in figure 4, which 

already contains the graph of discussed earlier; again the 

two cases are incompressible flow and compressible flow with an initial 
Mach number of 0.05. For 2-^1 and incompressible flow, the lateral 

component of turbulence is seen to increase steadily with 2^, in 

marked contrast to the decrease exhibited by the longitudinal component. 
The curve (of the lateral component) for compressible flow begins tq 
differ sensibly from the curve for incompressible flow for downstream 
Mach numbers above 0.3; above sonic speed compressibility is seen to 
effect a complete reversal of the curve. The over-all effect of com- 
pressibility on the contraction effect is thus much greater for the 
lateral than for the longitudinal component of the turbulence. 

The graph of equation (35) in figure 4 is primarily illustrative; 
a form more useful for engineering applications is given in figure 6. 

The single curve provides the variation of V v b 2 / v a 2 with both 2^ 
and c; 2]^ and e may be determined from the initial and final Mach 
numbers by means of equations (30) . 


DECAY CONSIDERATIONS 
Criterions for Negligibly Small Decay 

The basis of the present analysis of the contraction effect is 
embodied in equations (lla) relating the pre-contraction and post- 
contraction vorticity distributions. The simplicity of this result 
and its derivation arises from the neglect of the turbulent decay;’ by 
decay is meant the viscous dissipation and all the (nonlinear) inter- 
mixing processes of the eddies which together cause the mean turbulent 
intensity to diminish with time. The postulation of an inviscid fluid 
eliminated the viscous dissipation, and the limitation to very weak 
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turbulence eliminated the intermixing processes. (While there can be 
no dissipation in an inviscid fluid, the intermixing processes ordinarily 
associated with decay will occur.) In order to assess the influence of 
these assumptions equations (lla) will now be derived in a more general 
fashion with the Navier-Stokes equations as the starting point. For 
simplicity the fluid is taken to be incompressible, since the major con- 
clusions are unaffected thereby. 


General formulation of changes in vorticity . - By rearrangement and 
cross differentiation to eliminate the pressure term (reference 8, p. 578), 
the Navier-Stokes equations can be transformed into 


oo 

t-- 

to 

<u 


Do>]_ 

Dt“ 


dqi* 
“1 chc]~ 


dq.i* 
+ c °2 sr 


+ CO 


3q- 


'3 ctec-2 


2 

+ vy co-^ 


(36) 


and two similar equations, where co = , ooj is the vorticity and 

q* = q^’, q _2 ’, qj 1 is the resultant velocity. Now let q* be the sum 

of a stream velocity U, V, W and a turbulent velocity field 
£ = q^> 0-2 > < 3 . 3 > also, let curl U, V, W = 0, so that co is just curl q. 

Then equation (36) becomes (in tensor notation) 


Dco^ 

5b" 


m *t 

L_ 


SU ^3-1 

^ + “P ^ + 


_2 

vy 


(37) 


JL 


Contraction 


— — v — 

Decay 


and there are two similar equations. The first set of terms on the 
right-hand side is identified as the contraction effect, the second 
set as the decay effect. First the decay terms will be neglected in 
an attempt to recover equations (lla); then the neglected decay terms 
will be examined and criterions for their neglect arrived at. 

Neglect of the decay terms . - Equation (37) minus the decay terms 
reads, in expanded form, 


^(U+q-j_) 

Dt" = “1 ~£T- + 


ct(U+q 1 ) 

— 3 ^ + 

OX o *■ 


SOJ+q-L) 

5xZ 


(38) 


In this and the earlier equations — is the ’Lagrangian* derivative 
following the fluid motion. Now consider a line segment 5x^, Sxg, 8 x 3 
following the fluid motion: its Lagrangian derivative can be shown to be 
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D6Xj 

Dt 


= 5xi 


3(U+<li) 

OXi 


5x 2 


cEo 


+ Sx-3 


d(TJ-Hll) 

5E 


(39) 


and two similar equations. It can be seen that a solution of equa- 
tions (38) and (39 ) , together with their companion equations, is given by 

. ®1j ®2> “3 ~ 8x l> 8 x 2j 8x 3 

for all time t; this result is well known. Now complete the neglect of 
the decay terms by omitting the terms in q-^ in equation (38) and corres- 
pondingly in equation (39) . By this neglect the turbulent perturbations 
of the flow streamlines have been suppressed: this can be inferred from 

the revised equation (39). If the particles are at station A at a t im e 
t = 0 and reach station B at time t = t, there results 


and two similar equations. 


G>1. J 

*1 

But 


B Sx 


B 


Sx-^ 


5x 


B 


5*3 B 

5x. A 


5x-|_ A 


is just 


5 *2 E 

5xo A 


is l 


2> 


is l-z. Therefore equations (11a) have been recovered for the 


incompressible case (density ratio a = l). 

Consideration of inertial decay terms . - In equations (37) to (39) 
the decay terms not involving v are the inertial or intermixing terms. 
These are seen to be nonlinear. The condition for their neglect is 
evidently 


60, 


(40) 


and two similar conditions between q^ and V, q^ and W, respectively «, 
Tn a contraction like that of ,a wind tunnel the dominant velocity 

Su Bv cM 

gradients will be and these will be of the same order of 

absolute magnitude, A sufficient condition to replace (40) is therefore 




dU 

Sxp 

« 

^x 


(41) 


that is, all of the turbulent velocity gradients are very much less 
than the axial gradient of the stream velocity. This is essentially the 
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assumption underlying the distortion equations (10), which led directly 
to the vorticity changes (lla) in Taylor’s method. 


In statistical terms an approximate inference from equation (41) 
is, for isotropic turbulence, 



But by definition of X this may be written 



(42) 


The ’microscale’ X may be interpreted as a sort of average eddy diameter 
weighted in favor of the smaller eddies. Equation (42) may be' accepted 
as a practical criterion for the neglect of the inertial decay terms, 
equivalent to one of the two assumptions underlying equation (10) . The 
other assumption, neglect of viscosity, iB considered next. 


Consideration of viscous decay term . - The viscous decay term in 
equation (37) is the term containing v7 This term is linear and so will 
affect individual plane waves separately without mutual interference. 

The magnitude of the term may be estimated to a sufficient approximation 
by considerihg a wave carried along by the contracting stream and 
neglecting (for this term only) the distortion of the wave imposed by 
the contraction. Thus a component of the wave may be written 


“1 



k-jUt) 


Then, if the inertial decay terms of equation (37) are negligible, the 
equation reads, with 


cHJ 


dU 


dU 



dx 3 

« 

dx^ 


D^l dU 2 
dT = “1 35EI + V * “l 

= + v(k x 2 + k 2 2 + k 3 2 ) ©i 

( dU . J\ 

= w i y^ + v£ ) 
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Accordingly, viscosity may be neglected for. that portion of the spectrum 
which satisfies the inequality 


« 


dx^ 


(43) 


In the ‘initial period* of decay if the inertial decay criterion (42) 
is satisfied the major part of the spectrum will satisfy (43) . 


Rough Estimation of Mutual Effects 
of Decay and Contraction 

When the decay effects are not negligible compared with the c on - 
traction effects (see criterion developed in the last section) the 
theoretical basis of the present theory of the contraction effect is 
violated. Because negligible decay is more the exception than the 
rule, there is considerable incentive to attempt to apply the theory 
outside the valid range by means of assumptions concerning the simul- 
taneous effects of decay and contraction. 

Suppose, now, the decay and contraction are considered to occur 
alternatively in small steps, starting from isotropic turbulence. Each 
stream tube is considered to contract stepwise: between steps there is 

decay without contraction; at each step there is a sudden contraction 
without decay. Let the change in speed ratio per step be d2-]_, the 

reduction in u 2 due t o de cay be (du 2 ) D , and the reduction in u 2 
due to contraction be (du 2 ) c . Express the effect of decay in the 
absence of contraction in the form 



(44) 


where is a function of the time of travel (decay time) ' t, and 

the effect of contraction in the absence of decay in the form 



( 45 ) 


The corresponding differential forms are 
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D'Uq) 

^ir d2 i 



(46) 

(47) 


The assumption is now made that equation (46) applies to the decay effect 
per step and equation (47) to the contraction effect per step, the only 

o 

interaction being in the common u . The total effect per step is then 

— ' odi) + Whj 

u- u 



whence upon integration the over-all effect is 


u 1 2 * * * * 

= = cup DUp- (* 8 ) 

2 

U A 

That is, if the effect of contraction alone is expressed by C(Zp 
(equation (45)) and the effect of decay alone by D(Zp (equation (44)), 

then the joint effect under the assumption is stressed by the 
product C(Zp D(Zp. 

Equation (48) is intended to provide a very rough adjustment of 
the theoretical contraction effect C(zp to account for decay. This 


1 &\1 C p 

It is known that in the ’initial' period of decay ~ - u dt. 

_ u 2 ’ 

Equation (46) amounts to replacing the -u 2 on the right-hand side by 

-(u 2 )decay only’ 3ome defense may be made of this approximation, con- 

sidering the progressive deviation from, isotropy. 


2378 


NACA TN 2606 


35 


adjustment will be made in the attempt to compare the theory with experi- 
mental results in which the decay effects are of the same order as the 
contraction effects. 

Equation (48) refers to the longitudinal velocity component u; an 
equation of the same form is obtained for the lateral component v. For 
both cases the function D(2p) is taken to be the right-hand side of 

the empirical decay law (reference 9) 



for isotropic turbulence in the initial period. The decay time t(Zp) 

in the formula is the time required by a particle of the main stream to 
pass through the contraction, the initial velocity being U A and the 

final velocity 1]_U A . 


COMPARISON WITH EXPERIMENT 

There appear to have been no experimental investigations with 'which 
to compare the predicted changes imposed by a stream contraction on the 
spectrum of the turbulence, or on the correlation tensor of the tur- 
bulence. The available experimental data seem to be limited to measure- 
ments bearing on the changes in the root mean square velocity components. 
These data apply, moreover, to conditions outside the proper scope of 
the present theory in that large decay effects are present. The experi- 
mental data are therefore compared with a crude extension of the theory 
in which the decay is allowed for in first approximation. (See pre- 
ceding section.) 

The most extensive data are those of MacPhail, (reference 10) which 
in effect cover a range of contraction ratios from 2]_ = 1 to 2^ = 9*65 

inasmuch as measurements were made at various stations along the con- 
traction. Isolated points for particular contraction ratios were 
obtained from investigations made for other purposes by Dryden and 
Schubauer (reference 11 ) and by Hall (reference 12). Only those points 
were chosen for which the initial turbulence was indicated to be approx- 
imately isotropic. In the case of reference 11, data for the case of 
screens in the settling chamber were excluded because the final tur- 
bulence level was sensibly indistinguishable from the residual noise 
level. 
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In table I are listed, for the three experimental arrangements, the 
parameters used in the estimation of the decay factor (equations (48) 
and (49)). In reference 10 the Initial stream velocity and integral 

scale of turbulence were given. In reference 11 the value of 

was given, and the value of was taken to be 0.05 feet, the only 

scale mentioned; it was not clear, however, whether this value of scale 
applied with or without screens. In reference 12 the value of was 

inferred from collateral information and is somewhat uncertain; the 
scale was estimated from the dimensions of the honeycomb. In all 

three experimental arrangements the initial relative levels of turbulence 
were specified. The decay time t of the turbulence was computed as the 
time for a particle to traverse the contraction; the value arrived at 
for Hall’s data (reference 12) reflects the uncertainty in the 
assumed U^. 

Root mean square longitudinal velocity components . - The comparison 
of the theory, including estimated decay, with experiment for the longi- 
tudinal component of turbulence is given in figure 7. The theoretical 
curve, in each instance, is the product of a value computed for con- 
traction alone^ neglecting decay, (obtainable from fig. 5) and a 
value estimated for decay alone neglecting contraction. (See equa- 
tions (48) and (49).) The agreement with MacPhail’s data and with 
Hall's single point can be considered good. The agreement with the . 
Dryden-Schubauer point, on the other hand, is poor; a slight improve- 
ment would result on correction for the spurious contribution of the 
noise background . 

Root mean square lateral velocity components . - Comparison of the 
theory, again including estimated decay, with experiment for the 
lateral component of turbulence is given in figure 8. There is com- 
plete disagreement with MacPhail’s data and Hall’s single point, and 
on the other hand, good agreement with the Dryden-Schubauer single 
point. Thus there is the curious result that MacPhail’s and Hall’s 
data agree well with theory for the longitudinal component and disagree 
entirely for the lateral component, whereas the converse is true for 
the Dryden-Schubauer data. 

Discussion . - The uncertainty both in the manner of estimating the 
decay effect and in the data (table I) on which the estimate was based 
is still far from sufficient to account for the discrepancies between 
theory and experiment for the lateral component of turbulence . The 
very large amplification found by MacPhail is particularly hard to 
explain. On the other hand, the experimental data of the several 
observers show considerable disagreement, especially when differences 
in decay are allowed for. This disagreement would tend to cast doubt 
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on the validity of some of the data; the disagreement may also he in 
part a consequence, predicted hy the theory, of possible differences 
of the initial spectrums from each other and from isotropy. 


CONCLUDTOG REMAEKS 

The original aim of this paper was to provide a quantitative 
explanation of the observed changes in the root mean square velocity 
components of the turbulence of a wind-tunnel stream after passing 
through the tunnel contraction. The simplifying assumption of negligible 
decay was made to make the analysis tractable, although the decay and" 
contraction effects are ordinarily comparable. The analysis on this 
basis disclosed, in addition to the above integrated effects, pronounced 
changes in the spectrum of the turbulence. The changes in the shape of 
the spectral density curves, as distinguished from over-all changes in 
amplitude, would appear to be considerably less sensitive to modifica- 
tion by decay than would the mean square velocity components. For this 
reason, and because such spectral changes have not previously been dis- 
cussed, the emphasis of the present paper has been placed most heavily 
on these spectral effects. 

In particular, it has been found that the one-dimensional longi- 
tudinal spectrum, for isotropic turbulence exhibits a rather interesting 
change in shape downstream of the contraction; the center of gravity of 
the curve of spectral density versus longitudinal wave number is shifted 
substantially to higher wave numbers, the resulting distortion moving 
the peak of the curve well to the right of its initial position above 
the origin. The distortion is quite pronounced and would appear to be 
readily amenable to experimental observation. 

The restrictive assumption of negligible decay largely defeats, the 
original aim of the paper. Nevertheless, for practical reasons an 
attempt has been made to provide a crude extension to the theory in 
which decay is allowed for in first approximation. With this approxi- 
mation the theory^ has been compared with experimental values of the 
contraction effect on the longitudinal and lateral component root mean 
square velocity fluctuations. The agreement for the longitudinal 
component is good, whereas there appears to be almost complete disagree- 
ment for the lateral component, the experimental data themselves being 
in conflict. It is perhaps premature to attempt any general conclusion. 
For the present, the theory as augmented by the estimated decay effect 
may be useful in wind -tunnel-design applications. 

It is clear that the tunnel contraction effect on the components 
of turbulent intensity cannot be represented by fixed fractional changes 
independent of the character of the initial turbulence. Instead the 
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separate factors for the longitudinal and lateral components depend 
markedly on the spectrum of the turbulence . For initial isotropy, 
however, unique factors are predicted that, when decay is neglected, 
are independent of the details of the spectrum. 

Lewis Flight Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, August 30, 1951 
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APPENDIX A - SYMBOLS 


The following notation is used in this report: 

The subscripts, 1, 2, 3, refer to a rectangular coordinate system 
with the 1-axis alined with the axis of the main flow and directed down- 
stream, the 2-axis directed horizontally, and the 3-axis vertically,, 
Separate systems are used with the origins at stations A and B, respec- 
tively. (See fig. 1.) Vector and tensor notations are used inter- 
changeably; for example, k « k a = (k^/kg, k 3 ), where a = 1, 2, or 3, 

designates a vector with components k-^, kg, and k 3 . 

C(2^) function defined in equation (45) 

D (Z]_) function defined in equation (44) 

D edge length of cube within which turbulent field is defined 

e base of natural logarithms 

F a = F l> f 2 > or f 3 

F-j_ one -dimensional longitudinal spectral density (see equa- 

tion (7)) 

Fg , F 3 one -dimensional lateral spectral densities (see equation (7)) 


G(k) 

H 

Irn 
i = 

k 

k = 


function appearing in isotropic spectrum tensor 


constant 


(l k4oM 


Wo 

imaginary part of 

V-i 


dk 


K 3 = k^+k-L*, kg+kg*, k 3 +k 3 ‘, respectively 
amplitude of k -y/k 1 ^+kg^+k 3 2 j 
k a = (k]_, kg, k 3 ) wave number vector (station A) 


L 


longitudinal macroscale 
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V s 

*1 


M 

N 

Q = 

£ = 
Bop(£) 
Re 
r 

r = 

s 

t 

t 

U 

u,v,v = 


l l> l 2> 


stream velocity at station B divided by stream velocity at 
station A (see fig. l) 

stream breadth at station B divided by stream breadth at 
station A (see fig. l) 

stream height at station B divided by stream height at 
station A (see fig. l) 

Mach number of main stream 

amplitude of special isotropic spectrum tensor (see following 
equation (26)) 

Q a = (Q]_> %) disturbance wave amplitude vector 

q. a = (^l* *l2> £3) disturbance velocity vector 


correlation tensor (reference 6) 
real part of 

magnitude of r *= -y/r 1 2 +r 2 2 +r 3 2 

r a = ( r i> r 2> r 3) separation vector of two correlated points 

/ k 2 

parameter in equation (28) 1= _i_ + 

V T 2 



parameter in equation (28) 


(e-l)kl 2 


time 

main-stream velocity 

, Q.3 disturbance velocity components 


X length of wind tunnel contraction (distance between sta- 

tions A and B) 

used occasionally in place of Xj_ 


x 


2378 



NACA TN 2606 


41 


x = 
7 


e 

e o^7 

e 

X 

X = 



o 


'Z 

V 


<p 

0 ) = 
2 = 


x a = (Xp x 2 , x 3 ) position vector (station A) 
spectrum tensor (reference 3) 

constant in special isotropic spectrum tensor (see following 
equation (26)) (7 = l/L) 

contraction parameter (= X^ ^ l^\> see fig* l) 

alternating tensor defined after equation (9) 

polar angle (equation (31)) 

magnitude of x -J X^ 2 +X2 2 +X3 2 j 

x a = (x x , x^ } Xg) transformed wave number vector (station B) 
(’‘a = ^aAa) 

magnitude of 5 

£ a = (?i> £3) transformed position vector (station B) 

(see equation (10) and fig. 1(b)) 

summation over a for a = 1, 2 } 3 

stream density at station B divided by stream density at 
station A 

a volume 

viscosity 

azimuth angle (equation (31)) 
co a = (co-^ v> 2 , 0)3) vorticity vector 

2 = (2., 2 9 , 2*) vector amplitude of vorticity wave 

GG I Ci O 


Superscripts : 

A measured in vicinity of station A } upstream of contraction 

B measured in vicinity of station B, downstream of contraction 

* complex conjugate 
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Subscripts: 

A measured in vicinity of station A, upstream of contraction 

B measured in vicinity of station B, downstream of contraction 

take on values 1, 2, or 3 and designate tensor quantities 

1,2,3 specific values of a, p, 7 , or & 

A symbol with the mark ~ above it refers to a single plane wave. 

A bar over a symbol designates an average (usually a spatial average); 
a bar under a symbol designates a vector. 
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APPENDIX B 


EQUIVALENCE OF SPACE AND TIME AVERAGES IN STATISTICALLY STEADY, 

HOMOGENEOUS TURBULENCE 

The definitions of statistical homogeneity and statistical time 
independence -will first he made precise 0 Let F(x,y,z,t) he some 
property of a turbulent field that varies in time and from point to 
point; thus F may he the pressure, or any of the velocity components, 
or a correlation of velocity components at two points of fixed separa- 
tion, (x,y,z) being one of the two points. If, for all choices of the 
property F, (a) the average of F over a time T -* 00 is independent of 
(x,y, z), the turbulence is defined to he statistically homogeneous; if 
(b) the average of F over a volume V -*• 00 Is independent of t, the 
turbulence is defined, in the sense used herein, to be statistically 
steady or time-independent. The respective averages are supposed to be 
approached uniformly, in the mathematical sense, as T or V, 
respectively, approach infinity. (A statistically steady or "stationary" 
condition is defined differently in the theory' of random processes.) 

It will now be proved that if the turbulence satisfies the two / 
conditions (a) and (b), the time and space averages defined therein are 
equal. In this proof no resort will be made to the "ergodic hypothesis" 
of statistical mechanics, which leads to the equivalence of the time 
average and the "ensemble" average. The possibility of the joint 
existence of the conditions (a) and (b) probably amounts, however, to 
just as fundamental an assumption. 

The space average will be made over a parallelepiped of edges a, 
b, c and the time average over a time T, and then a limiting process 
will be applied. The average of F over both space and time is thus 



/ 


^b ( 

F . = Lim m 1 

1 



s,t a,t,c,T~ aboT 

J 

io ^ 



0 


0 v. 


F dx dy dz dt 


(Bl) 


Any order of integration is permissible, since the integration 
limits are constants. If the time integration is performed first the 
expression may be written 


F + = Lim 1 ' 

• 8,t a >' b ; c ; T-> ® abcT 

^c j 

^b 

-Ny 

~ T \ 
F dt ) 


0 

6o 

o\u 

0 J 
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By virtue of the postulated uniform convergence of the time and space 
averages the operation may he brought under the integral sign: 


Lim 


1 

* abc 


'"'c 

I 

v 

1 Lim i 

I 

\ 1 
F dt dz d y dz = 

r '.c 

"'b 

0 c 

0 ^ 

k \ 

0 J J 




(B2) 

dx dy dz 


where F^ is the time average of F. But, by condition (a), F t is 
independent of x ,. y, and z. Therefore 


F s,t = F t 


(B3) 


Alternatively, the space integration and limiting process may be 
performed first: 


s = Llm 1 
s,t T-** T 



Lim _1_ 
a,b,c r*» abc 



“>b | 


r* rj -1 



\ Lim 1 




F dx dy dz dt - 


- 0 v > 

0 ^ 

0 / ^ 

0 


F s dt 


(B4) 


where F a is the space average of F. By condition (b), F s is " 
independent of t; therefore, 


F j. = F 
s,t r a 

Equation (B3) and (B5) together state that 


(B5) 


F s = F t - F s,t (B6) 

or the space average, the time average and the space-time average are 
all equal. 
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TABLE I - BATA POE ESTIMATION OF BECAY 



n 

U A 

(ft/sec) 

L A 

(ft) 

t 

(sec) 

Vs? 

(ft/sec) 

MacPbail (reference 10) 

1.20 

3.55 

0.012 

0.19 

0.149 


1.60 

3.55 

.012 

.34 



2.55 

3.55 

.012 

.45 



4.90 

3.55 

.012 

.51 



9.65 

3.55 

.012 

.57 


Bryden-Sch.ubau.er 

6.6 

6.86 

.05 

1.31 

.114 

(reference 11 ) 






Hall (reference 12) 

5.2 

1.54 

.025 

1.22 

.046 
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Figure 5. - Variation of relative root mean square longitudinal velocity fluctuation with both speed ratio 

for axisymmetrlc contraction and isotropic turbulence at U - 1 . 
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Figure 7. - Comparison of predicted axisymmetric contraction effect with experiment for longi- 
tudinal component of turbulence, with decay allowed for in first approximation. Initial' 
isotropic turbulence assumed. 
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Figure 8. - Comparison of predicted axisymmetric contraction effect with exper 
component of turbulence, with decay allowed for in first approximation. Ini 
turbulence assumed. 



